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Resumen 

En el estudio de poblaciones de mamiferos marinos, es comun el uso de censos visuales 
para la clasificacion de las zonas de distribucion de una especie; la secuencia de los resultados 
de estos censos en el tiempo, puede ser considerada como un proceso estocastico. A1 analizar 
estos resultados en conjunto con otras variables mediante analisis multivariados, su ordenamiento 
genera hipotesis donde los avistamientos son variables causales y/o de respuesta a n varia¬ 
bles ecologicas, consecuentemente estas hipotesis pueden ser validadas por modelos estocas¬ 
ticos. El modelo cadenas de Markov se aplico para validar la hipotesis la cual postula que la 
distribucion espacial del manati del Caribe depende de las variaciones drasticas en la intensi- 
dad del viento. El modelo asigno una probabilidad de 0.40 de que aumente el numero de ma¬ 
nages despues de un cambio drastico en la intensidad del viento, de 0.56 de que disminuya y 
de 0.04 de que no vane. La validez del modelo dependera del comportamiento de las variables 
ecologicas analizadas en el tiempo, siendo valido si las variables analizadas tienen un com¬ 
portamiento semejante entre anos; por lo tanto, la distribucion de probabilidades tendra un 
valor de significacion p igual al valor de analisis del comportamiento de las variables y el es- 
tadistico asociado al valor p correspondera al de la prueba de contestation. 


Abstract 

In the population study of marine mammals, it’s common the censuses visual use for the 
zones classification of species distribution; the results sequence of these censuses in the time, 
could be considered like a stocastic process. Upon analyzing these results on the whole with 
other variables by means of multivariate analysis, their ordination generates hypothesis 
where they are causal variables and/or answer to n ecological variables, consequently these hy¬ 
pothesis could be validated for stocastic models. The chains Markov model worked in order 
to validate the hypothesis which it postulate that the spatial distribution of the Caribbean 
mannate relies on the drastic variations in the wind intensity. The model assigns a probability 
of 0.40 that it increase the mannate number after drastic change in the intensity of the wind, 
of 0.56 that it diminish and 0.04 of that it don’t vary. The validity of the model will rely on the 
behavior of the ecological variables analyzed in the time, being been worth if the analyzed va¬ 
riables have a similar behaviour between years; therefore, the probability distribution a sig¬ 
nificance value of p same to the value of the variables behaviour analysis and the associate sta¬ 
tistic to the p value will correspond to the contrastation test. 


Palabras claves: Cadenas de Markov, Procesos estocasticos, Manati, Caribe Mexicano. 
Key words: Markov’s chain, stocastic processes, mannate, mexican Caribbean. 


11 


Axis-Arroyo, J. et al. 


INTRODUCTION 

En ecologfa de comunidades los procesos sinecologicos, aun en las llamadas comu- 
nidades climax, son procesos con un nivel de incertidumbre asociado a ellos, pero que 
pueden presentar un comportamiento en espacio y tiempo con un cierto grado de regula- 
ridad (Gauch, 1982; Jongman et al., 1987), esto permite describir a las variables gene- 
radoras del proceso mediante un modelo probabilfstico que las considere en forma cuan- 
titativa y no cualitativamente. Cuando las variables que determinan un proceso sinecolo- 
gico se observan en espacio y/o tiempo fijo y evolucionan en esas mismas dimensiones, 
el proceso sinecologico es estocastico (Legendre y Legendre, 1984). 

En el estudio de poblaciones de mamfferos marinos, es comun el uso de censos vi- 
suales o de algun fndice de abundancia basado en estos para la clasificacion de areas. La 
presencia de individuos de una especie determinada en un espacio dado en el tiempo, 
puede ser descrita facilmente por un modelo probabilfstico con variaciones temporales. En 
este contexto, una secuencia de avistamientos en el tiempo puede ser considerada como 
un proceso estocastico; si los datos de avistamiento se analizan conjuntamente con otras 
variables mediante cualquier analisis multivariado, el ordenamiento de estas variables 
generarfan hipotesis sobre procesos en los que los avistamientos seran la variable causal 
y/o de respuesta, en consecuencia estas hipotesis podrfan ser validadas utilizando mode- 
los estocasticos. 

En 1996 se analizaron las variables que determinan el uso de habitat del manatf del Ca- 
ribe (Trichechus manatus manatus) en la costa oeste de la Bahia de Chetumal, Quintana 
Roo, utilizando el analisis de direccion (Axis-Arroyo, 1998). Este metodo multiva¬ 
riado, evalua mediante un modelo generalizado el total de relaciones causales que describen 
un sistema, pero no las relaciones que se presentan entre las variables. Considerando 
solo el ajuste de los valores de estas relaciones con la estructura de datos propuesta por el 
modelo, los valores obtenidos de las relaciones entre variables no tienen robustez estadfstica 
por si solos, pero pueden generar hipotesis sujetas a comprobacion. 

En el presente trabajo se aplicara el modelo estocastico de cadenas de Markov para eva- 
luar una hipotesis generada por el analisis multivariado, en el que los avistamientos de una 
especie en un area son la variable de respuesta a n variables ecologicas. 


MATERIALES Y METODOS 

Al analizar las relaciones entre las variables, una de las hipotesis obtenidas, fue que la 
presencia de un manatf en areas de aguas abiertas (sin resguardo), en dos de las tres tem- 
poradas climaticas en la zona (nortes y secas), depende de las variaciones drasticas de la 
intensidad del viento. Esta hipotesis sera la que evaluaremos mediante el modelo esto¬ 
castico de cadenas de Markov, para lo cual utilizaremos los datos de los muestreos reali- 
zados en los dfas con una intensidad del viento 25% mayor a la intensidad promedio del 
mes, comparandolos con el muestreo inmediatamente anterior realizado en dfas sin va¬ 
riaciones drasticas, durante los meses de octubre del 94 a junio del 95; en la Tabla 1 se re- 
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portan los datos del numero de manatfes avistados en el total de las areas muestreadas y 
que serviran como base para el presente analisis. Para el calculo de las probabilidades de 
transicion los avistamientos entre muestreos se compararon por cada area y no en total. 

Un proceso estocastico es un conjunto referido de variables aleatorias {Xt} en donde 
el subfndice t toma valores de un conjunto T y Xt representa una caractenstica medible en 
el tiempo t (Gnanadesikan, 1977). 

Un proceso markoviano presenta las siguientes caracteristicas: 

a) Es un proceso discreto en el tiempo (t = 0, 1,2,..., n). En nuestro caso, esta definido 
por cada uno de los muestreos. 

b) Es un proceso con un numero finito de estados mutuamente excluyentes. Consi- 
derando que el parametro de analisis es el avistamiento de animates, el proceso tendra 3 
posibles estados en relacion al tiempo anterior: 

- Mismo numero de avistamientos (estado estacionario Pi). 

- Aumento en el numero de avistamientos (Pj). 

- Disminucion en el numero de avistamientos (Pk). 

c) Esta gobernado por las probabilidades de transicion Pij (el paso del estado i al j 
esta dado por la probabilidad Pij). La transicion de estados se determinara por la proba- 
bilidad de que se mantenga, aumente o disminuya el numero de avistamientos entre los tiem- 
pos t y t+1. 

d) El estado j del proceso depende solamente del estado i y es independiente de la 
historia del proceso. El estado en el tiempo i dependera exclusivamente del numero de avis¬ 
tamientos registrados en el tiempo i-1. 


Muestreo 

Relacion 

Muestreo 

Relacion 



entre 



entre 

Anterior 

Posterior 

muestreos 

Anterior 

Posterior 

muestreos 

5 

3 

Disminucion 

4 

3 

Disminucion 

5 

12 

Aumento 

8 

1 

Disminucion 

2 

6 

Aumento 

2 

1 

Disminucion 

4 

2 

Disminucion 

4 

6 

Aumento 

7 

6 

Disminucion 

10 

6 

Disminucion 

7 

12 

Aumento 

5 

1 

Disminucion 

4 

1 

Disminucion 

6 

8 

Aumento 

8 

3 

Disminucion 

1 

2 

Aumento 

7 

7 

Sin cambio 

1 

3 

Aumento 

9 

4 

Disminucion 

2 

3 

Aumento 

5 

3 

Disminucion 

4 

1 

Disminucion 


Tabla 1 . Numero de manatfes avistados en el total de las areas muestrales, registrados antes y 
despues de variaciones drasticas en la velocidad del viento (>25% de la velocidad promedio de 

los tres dfas anteriores). 

Table 1. Mannati number watched in the all areas show them , registered before and after drastic 
variations in the wind speed (>25% average speed of the three previous days). 
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e) Tiene un conjunto de probabilidades iniciales P {X0=i} para toda i. El conjunto 
de probabilidades iniciales de determinara por la diferencia entre el primer y segundo 
muestreo. 

Al plantear la evaluacion de la hipotesis como un proceso Markoviano es conve- 
niente especificar los dos conceptos basicos del comportamiento de un sistema en una 
cadena de Markov: 

a) Estado del sistema: Los valores que describen por completo la position del sis¬ 
tema en cualquier instante del tiempo son importantes para el comportamiento futuro del 
sistema; en este caso, el numero de avistamientos registrados en cada una de las areas en 
cada muestreo; b) La transicion de estados (Protter y Morrey, 1969). Al tratarse de 
tiempos discretos la transicion de estados son las variaciones en el numero de avista¬ 
mientos registrados entre los tiempos t y t , lo cual representara la evolution del sistema 
en el tiempo. 

A cada posible transicion del estado i al j se asocia una probabilidad P... que en pro- 
cesos Markovianos se define como la probabilidad de transicion de un paso, con valores 
0<P..<1 (tomara valor 0 cuando no sea posible realizar una transicion del estado i al j, y 
de 1 si estando en el estado i pasa exclusivamente al estado j), lo cual nos permite definir 
probabilidades de transicion a partir de un estado i a un estado j. Si tenemos n estados ob- 
tendriamos una matriz P denominada matriz de transicion, donde: 

P P — P ~~ 

i i i un 

■ 

p p —p 

. m in nun _ 

Considerando los 3 estados del presente analisis, nuestra matriz de transicion estarfa 
dada por: 

P = 


P 

u 

P. 

K 


P 

P' 

P 

kj 



Donde: 

P. = No hay cambio en el numero de avistamientos despues de un estado estacionario. 
P..= Aumento en el numero de avistamientos despues de un estado estacionario. 

P. k = Disminucion en el numero de avistamientos despues de un estado estacionario. 
P..= Aumento en el numero de avistamientos despues de un aumento. 

P.= No hay cambio en el numero de avistamientos despues de un aumento. 

P. k = Disminucion en el numero de avistamientos despues de un aumento. 

P Rk = Disminucion en el numero de avistamientos despues de una disminucion. 

P ki = No hay cambio en el numero de avistamientos despues de una disminucion. 

P k = Aumento en el numero de avistamientos despues de una disminucion. 
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Esta matriz de transicion cumple con las siguientes caracterfsticas del proceso: 

a) En cualquier momento cualquier manati puede estar en cualquier area; por lo que en 
un area x puede observarse cualquiera de los tres estados, es decir los estados son proba- 
bilfsticos. 

b) La probabilidad de que aumente, disminuya o se mantenga el numero de manaties 
avistados en un area; es relativa al numero de avistamientos en el tiempo inmediata- 
mente anterior, es decir cumple con la propiedad Markoviana. 

c) La suma de las probabilidades de los 3 eventos de cada renglon de la matriz de 
transicion sera 1, es decir tiene a 1 como eigenvalor y ninguno de ellos excede a 1 en va¬ 
lor absoluto. 

Considerando los avistamientos entre muestreos por cada una de las areas, la diferencia 
entre el primer y segundo muestreo no tendra un antecesor, por lo que se considerara como 
el primer eslabon de la cadena markoviana, al no tener antecesor la matriz de transicion 
tendra solo 3 probabilidades: 

p =i p , p j P J 

La matriz sera un vector al que se le denominara vector inicial de distribution {X () }, 
con componentes no negativos y sumatoria 1 y que tomara los valores de la siguiente 
matriz: 


P. = [0.02 0.27 0.71] 

Las diferencias entre los estados 2° y 3° y subsecuentes nos generaria la matriz de 
transicion asociada al proceso y que en nuestro caso seria: 


0.00 

0.37 

0.63 

0.01 

0.31 

0.68 

0.06 

0.43 

0.51 


Por lo tanto, esta matriz modela las probabilidades de que se observe cualquiera de las 
9 transiciones generadas por los 3 estados. Considerando el tiempo, estas probabilidades 
de transicion deben ser estacionarias para poder deducir la probabilidad de que el sis- 
tema se encuentre en el estado j despues de n transiciones, a partir del estado inicial i, lo 
cual se denota por P.. (n) , lo cual es una probabilidad de transicion de n pasos. 

Para que el analisis del proceso Markoviano tenga sentido en el tiempo la matriz de tran¬ 
sicion debera de cumplir con las siguientes propiedades: 

1) Ser irreductible: que sea posible pasar del estado i al j en algun numero finito de tran¬ 
siciones para los estados i y j, esto es: 

P.. (n) >0 para algun (n) para toda ij 
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2) Ser aperiodica: un estado i tiene un periodo t>l si es posible para la cadena estar en 
el estado i unicamente en multiplos de t, es decir: 

P (n) = 0; siempre que n no sea un multiplo de t 

Si una cadena markoviana es finita, irreductible y aperiodica es una matriz ergodica y 
en consecuencia regular (Morrison, 1976). 

Cuando una cadena de Markov es regular la distribucion probabilfstica tiene a una 
distribucion limite comun conforme n aumenta, independientemente del estado inicial 
en que empezo el proceso y se determina por: 

limite ^ P.. (n) = p.; j = 1, 2, n 

n®a ij *j ° ’ 

Donde p. satisface las ecuaciones del estado estable: 
j 


n >o 

N 



N 


n=lR P.. ;j = 1,2,N 

j = l 

Esta ultima formula expresada en notation matricial toma la siguiente forma: 

n = np ,) 


La probabilidad p. es la probabilidad del estado estable y representa las probabilida- 
des de que el sistema se encuentre en cada estado al principio o fin de una transicion, 
despues de que haya ocurrido un numero lo suficientemente grande de ellas como para lo- 
grar una estacionalidad en las probabilidades. 

Considerando que: 

La n-esima potencia de la matriz P se denota como: 

P ={P.. (n) } 

n ij J 

El vector de distribucion en el inicio del proceso es dado por: 

X <n) = [X ( (n), X 2 (n),X(n)] 
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Y el vector de distribucion en el inicio del proceso es dado por: 

X (0) = [X l (0), x 2 (0), X n (0)] 

La relacion entre los vectores de distribucion en el n-esimo periodo y en el inicio del 
proceso se estableceria por la formula siguiente: 


X< n > = x * 0) P n 


En consecuencia, la distribucion de probabilidades en el estado limite es determinada 
por: 


Xf a) =Xf 0) L 2) 


Donde: 

X (a) = Distribucion de probabilidades del estado limite. 

X (0) = Distribucion de probabilidades en el inicio del proceso. 

L = Matriz de transition limite. 

En nuestro ejemplo X (0) toma los valores del vector P. y la matriz L los de la matriz P r 

RESULTADOS 

Los estados y posibles transiciones del sistema se especifican en el siguiente diagrama: 



En donde: 

Los circulos representan los estados. 
i = Mismo numero de avistamientos. 
j = Aumento en el numero de avistamientos. 
k = Diminution en el numero de avistamientos. 

A1 aplicar la formula 2 se obtuvo la distribucion limite de probabilidades: 

X< a) = [0.04 0.04 0.56] 


Lo cual representa que: 

Hay una probabilidad 0.40 de que aumente el numero de manaties despues de un au¬ 
mento drastico en la intensidad del viento. 

Hay una probabilidad de 0.56 de que disminuya el numero de manaties despues de 
un aumento drastico en la instensidad del viento. 
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Hay una probabilidad de 0.04 de que no vane el numero de manaties despues de un au- 
mento drastico en la intensidad del viento. 

DISCUSION 

El supuesto mas cntico del proceso markoviano en la presente aplicacion (y considera- 
mos que en todas las aplicaciones en ecologia), es la condicion de la estacionalidad de las pro- 
babilidades de transicion; esto es, que su matriz P permanezca igual en todo el periodo. 

Al considerar tal condicion, la validez del modelo dependera de las diferencias de las 
variables ecologicas en el tiempo, que determinaron el proceso estocastico. Como en 
este ensayo se consideraron las temporadas de Nortes y secas, el modelo describira las po- 
sibilidades de transicion en tales temporadas y tendra validez en el tiempo si las variables 
analizadas tienen un comportamiento estable entre largos periodos. 

Estadisticamente el modelo sera valido si no existen variaciones significativas entre el 
comportamiento del total de las variables en los tiempos t y t+1; por lo tanto, la distribu- 
cion de probabilidades tendra un valor de significancia p igual al valor del analisis del 
comportamiento de las variables y el estadfstico asociado al valor p corresponded al de 
la prueba que contrasta el comportamiento de las variables (Linares et al., 1986). 

Si el comportamiento de las variables no tienen diferencias estadisticamente signifi¬ 
cativas en el tiempo, las distribuciones de probabilidad seran las mismas. Si se presentan 
diferencias, estas seran un indicador de perturbaciones en la dinamica del sistema y de esta 
manera, el numero de avistamientos tendra un valor parametrico que cuantifique las per¬ 
turbaciones que ocurran en el medio. 

Considerando las variables en forma individual, un estado tiene un valor probabilistic 
generado por un conjunto de datos en los que tIT; por lo que las variaciones de cada t al- 
teran el proceso estocastico, la estabilidad de las condicionantes causales de estado de- 
penderan del apego de sus distribuciones a la distribution normal; sin embargo, la 
misma naturaleza estocastica del modelo da un intervalo de tolerancia, por lo que consi- 
deramos que una condicionante causal de estado es confiable si los datos de su distribu¬ 
cion cumplen con el teorema de Tchebichev. 
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